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Abstract. We consider the trajectories of points on under sequences of 
certain folding maps associated with reflections. The main result character¬ 
izes collections of folding maps that produce dense trajectories. The minimal 
number of maps in such a collection is d -|- 1. 


1. Introduction 

The subject of this paper is the reconstruction of fuii radiai symmetry from 
partiai information. A function on is radiai if and oniy if it is symmetric under 
reflection about arbitrary hyperpianes through the origin. Since the orthogonai 
group 0{d) acts transitiveiy and faithfuiiy on an equivaient statement is that 
the reflections generate 0{d). 

It is well-known that a finite set of reflections suffices to generate a dense subgroup 
of 0{d). In the plane, the composition of two reflections is a rotation by twice the 
enclosed angle. If the angle is incommensurable with tt, then the multiples of the 
rotation are dense in the circle. Likewise in dimension d > 2, the reflections at d 
hyperpianes in general position generate a dense subgroup of the orthogonal group. 
From any given starting point, the subset of points that can be reached by a suitable 
composition of these reflections is dense in the centered sphere that contains the 
point. The corresponding random walk is almost surely equidistributed on the 
sphere. 

Here we study the corresponding issues for a family of piecewise isometries that 
fold the sphere across a hyperplane onto a hemisphere. We label a folding map by 
the unit vector in the direction of the target hemisphere, and consider sequences 
of such maps indexed by a given set of directions G C The principal aim of 

this article is to answer the question: When can a dense subset of points in 8"^“^ be 
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reached from an arbitrary starting position by applying folding maps with directions 
chosen from G? The question came up in prior work on the convergence of random 
sequences of symmetrizations [1]. Our motivation will be discussed in Section 3. 

It turns out that a pair of obvious necessary conditions on the set of directions 
— one geometric and one algebraic — is also sufficient (Theorem 1). If a set of 
directions meets the conditions, then the only functions on that increase under 
composition with each of the corresponding folding maps are radial. Furthermore, 
the random walk generated by randomly alternating these maps is uniquely ergodic 
(Theorem 2). The invariant measure gives positive mass to all non-empty open 
subsets, but does not agree, in general, with the uniform measure on 


2. Main results 

We begin with some definitions. Let be the standard sphere, viewed as 
the set of unit vectors in The geodesic distance d{x,y) on given by the 

enclosed angle between x and y, is related to the chordal distance in by 

\x — y\^ = 2 — 2 cos d{x, y). 

The symbol 0{d) refers to the group of orthogonal linear transformations on R'^, 
and SO[d) to the orientation-preserving subgroup. 

The uniform probability measure on the sphere induced by Lebesgue measure on 
R"^ is denoted by tr. It will be used as a reference measure throughout the paper. By 
a null set we mean a subset A C with o'(A) = 0. All sets and functions under 
consideration are understood to be Borel measurable, and measures are assumed to 
be regular Borel measures. 

For a direction u G let = x — {2x • it) u be the reflection of a point 

X G R"^ about the orthogonal hyperplane if*-. Clearly, Ru = R-u- Let Hu = {x G 
R"^ I X • It > 0} be the open positive half-space associated with it, and define 


Fu{x) 


X , if X G Hu , 
RuX , otherwise. 


We call Fu the folding map in the direction of it, as it folds M'* across the crease 
at If*- onto the closed nonnegative half-space. It is nonlinear and idempotent = 
Fu). Since is isometric, Fu is non-expansive 


|F„(x) - Fu{y)\ < |x - 1/1, 

with strict inequality if x and y lie on opposite sides of it-*-. 

The two-point symmetrization of a real-valued function (f on is the 

equimeasurable rearrangement defined by 

|max{(()(i?„x),(()(x)}, AxGHu, 

bu(P\x) = < 

1 min{(()(i?„x), (()(x)} , otherwise. 

Note that Suif = if a-nd only A (j) o Fu > 4>. 

Let (un)n>i in be a sequence of directions in We study the trajectory 
of a point x in or 8"^“^ under the corresponding sequence of folding maps, given 

by 


( 1 ) 


Xo = X 


Xn = Fu^Xn -1 forn > 1. 
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Under what conditions on the sequence (m„) are these trajectories dense in 
Our main result shows that it suffices to choose the directions from a small subset 

Theorem 1 (Transitivity). Let G C be a set of direetions. If 

(Cl) the open half-spaces Hu with u G G cover , and 

(C2) the reflections {Ru \ u € G} generate a dense subgroup of 0{d), 

then there exists a sequence (u„)„>i in G such that for every starting point x G 
the trajectory defined by Eg. (1) is dense in 

Both assumptions on G are clearly necessary for the existence of even one dense 
trajectory. In Section 4, we show that the geometric condition (Cl) is equivalent 
to the origin lying in the interior of the convex hull of G in (Proposition 4.1(c)). 
In particular, G must contain at least d + 1 directions that span Condition 
(C2) can be replaced by an explicit sufficient condition (Proposition 4.2). Any set 
of directions that spans R'^ can be augmented by one more direction to satisfy both 
conditions (Proposition 4.3). 

As for the conclusion, we claim that there is a sequence of directions that produces 
dense trajectories simultaneously for all starting points on the sphere. A useful 
consequence of Theorem 1 is the following characterization of radial functions in 
terms of two-point symmetrizations. 

Corollary 1. Let (j> be a continuous function on R'^, and let G C 8“^“^ satisfy (Cl) 
and (C2). Then 


Su4> = 4> for all u G G (f is radial. 

The corollary improves upon the statement that for any non-radial function there 
exist two-point symmetrizations that push the function towards its symmetric de¬ 
creasing rearrangement. This observation appears as a key step in Baernstein and 
Taylor’s proof of Riesz’ rearrangement inequality on the sphere [2, Theorem 2 and 
p. 252], and many other classical results (see also [3, p. 226], [4, Lemma 6.5], [5, 
Lemma 2.8]), and [6, Lemma 3.11]. The corollary extends directly to measurable 
functions modulo null sets. 

For the proof of Theorem 1 in Section 5, we consider subsets of the sphere that 
are positively invariant under the folding maps indexed by a set G. Assuming G 
satisfies (Cl), every positively invariant subset of the sphere is invariant, up to null 
sets, under the corresponding reflections (Lemma 5.1). Condition (C2) precludes 
the existence of non-trivial invariant compact subsets of the sphere (Lemmas 5.2 
and 5.3). For any given starting point, the union of all possible trajectories under 
sequences of folding maps indexed by G is a positively invariant subset. We argue 
that this union is dense, and then select the desired dense trajectory from it. Our 
construction yields a sequence of directions (u„) with the following property: For 
any e > 0, there exists a number N such that every trajectory intersects every ball 
of radius e in within the first N steps, see Eq. (5). 

In the last two sections, we consider sequences of directions {Un)n>i that are 
chosen independently at random from a probability distribution p, on that is, 

N 

P{U„ G AnfloT n=l,...,N) = /i(A„) 

n—1 
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for any Borel sets Ai,..., C The trajectory defined by 

(2) Xq = X, Xn = Fu„Xn-i for n > 1 

will be called the random walk generated by ^ with starting point x. We are most 
interested in examples where /r is concentrated on a finite set of directions. 

Let G C be the support of p,, that is, the smallest closed subset of full /r- 
measure. If G satisfies the hypotheses of Theorem 1, then the random walk is almost 
surely dense (Corollary 5.5). Our second theorem strengthens this observation. To 
state the result, define a linear transformation on measures by 

(3) T^^v{A)= [ ,y{F-^{A))dfi{u) 

for A C and every Borel measure ly. By definition, maps the probability 

distribution of to the distribution of Xn+i for each n > 0. 

Theorem 2 (Unique ergodicity). Let ^ be a regular Borel probability measure on 
whose support satisfies (Cl) and (C2), and let be given by Eg. (3). There 
is a unique regular Borel probability measure p on with = p. The support 

of p is the entire sphere . 

We refer to p as the invariant measure for the random walk. By uniqueness, 
it assigns zero or one to every invariant subset. The invariant measure governs the 
behavior of the random trajectories {Xn) through Birkhoff’s ergodic theorem. 

Corollary 2. Let A C be a Borel set. Under the assumptions of Theorem 2, 
lim ^ \ {n < N \ Xn G A} \ = p{A) 

A—>oo JV 

almost surely for p-almost every starting point x S 

Thus, typical trajectories are equidistributed according to the invariant mea¬ 
sure p. Since p gives positive measure to every non-empty open subset, this provides 
another proof that almost all trajectories are dense. 

The proof of Theorem 2 is given in Section 6. To construct the invariant measure, 
we study the Markov chain associated with the random walk from Eq. (2) through 
its adjoint action on functions, defined by 

finix) = E^{(j){Xn)) , 

where x is the starting point of the random walk. Under the assumptions of The¬ 
orem 2, we show that for each continuous function (j) on the sphere, the sequence 
(pn converges uniformly to a constant 0 (Proposition 6.2). Writing 0 = J (p dp for a 
suitable measure p on the sphere, we show that p is uniquely determined by the map 
0 I—>■ 0, and invariant under the transformation Since 0 lies strictly between 

the maximum and minimum of 0 unless 0 is constant, we find that p{A) > 0 for 
every non-empty open set A. 

In Section 7, we consider the relationship between the invariant measure p and 
the uniform probability measure a on If the directions {Un) are uniformly 

distributed on the sphere, then p is uniform as well. In that case, Corollary 2 implies 
that the random trajectories are almost surely equidistributed on the sphere. More 
generally, the uniform measure is invariant if and only if the distribution of Un is 
even under the antipodal map u i—>■ —u (Proposition 7.1). It is an open question 
whether the invariant measure is always mutually absolutely continuous with respect 
to the uniform measure. 
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3. Motivation and related work 

Symmetrizations are equimeasurable rearrangements of sets and functions which 
are commonly used for proving that certain optimization problems have radially 
symmetric solutions. For instance, the sharp constants in the Young, Sobolev, 
and Hardy-Littlewood-Sobolev inequalities are assumed among radial functions, 
the ground state of the hydrogen atom is symmetric decreasing, and the membrane 
of given area with the lowest fundamental frequency is shaped like a disk. 

A standard technique for establishing geometric inequalities is to approximate 
full radial symmetrization by a sequence of simpler symmetrizations. Two-point 
symmetrization has been used in this way to prove the isoperimetric inequality on 
spheres [8], and sharp inequalities for path integrals [5, 6, 7]. The convergence of 
random sequences of symmetrizations has received some attention in the literature, 
most notably in the work of Klartag [9] on Steiner symmetrizations of convex sets. 
Convergence of two-point symmetrizations is less well studied. Very recently, De 
Keyser and Van Schaftingen have considered random symmetrization processes with 
time correlations [10]. Open questions in the area include precise conditions for 
convergence, and bounds on the rate of convergence. 

To explain the relationship with our results, let </> be a nonnegative continu¬ 
ous function with compact support on Consider a sequence {4>n) of functions 
obtained from tp by iterated Steiner symmetrization according to a sequence of di¬ 
rections. Any limit point ip of {pn) has at least some reflection symmetries, that is. 
Ip o = Ip for all u in a non-empty set of directions G C If G, which is de¬ 

termined by the sequence of Steiner symmetrizations, satisfies (C2), it follows that 
Ip is radial. Then one can identify ip as the symmetric decreasing rearrangement 
of (p and conclude that the entire sequence of symmetrizations converges to ip [1, 
Corollary 2.3b]. This provides a sharp condition for an i.i.d. sequence of Steiner 
symmetrizations to converge to the symmetric decreasing rearrangement. 

If, instead, pn is obtained by iterated two-point symmetrization at hyperplanes 
in R'^ that do not contain the origin, then any limit point ip of the sequence has the 
property that Sup) = P, that is, ip o > ip for all u in a some non-empty set of 
directions G C If G satisfies (Cl) and (C2), then according to Corollary 1, 

the function ip is radial. In combination with [1, Theorem 2.2] this yields a sharp 
condition for the convergence of i.i.d. sequences of two-point symmetrizations, 
the first result of this type. Our hope is that a full analysis of the random walk 
defined in Eq. (2) would shed light on the best achievable rate of convergence. We 
note that for non-convex sets, the strongest known lower bound on the rate of 
convergence of Steiner symmetrizations was proved by comparison with two-point 
symmetrization [1, Corollary 5.4]. 

There are many known results analogous to Theorem 1 for other collections of 
maps associated with linear isometries of spheres. For example. Crouch and Silva 
Leite [11] found pairs of one-parameter subgroups of SO{d) which generate all of 
SO{d), and produced upper bounds on the number of elements from each subgroup 
required to generate any element (see also Levitt and Sussmann [12] for related 
results). An example of this phenomenon is the Euler angles decomposition, a 
formula for writing any element of SO{3) as a product of three rotations about the 
X and y axes. Rosenthal [13] has established bounds on the rate of convergence to 
the steady state for random walks generated by conjugacy classes of planar rotations 
on SO{d). Porod [14] has obtained analogous results for random walks on 0{d), 
U{d) and Sp{d) generated by the conjugacy class of reflections. 
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All of the above results can be translated into statements about trajectories on 
via the standard actions of the groups. Historically, they were preceded by 
results on the symmetric group Sm with transpositions playing the role of reflec¬ 
tions. Dixon [15] proved that the probability that two randomly selected elements 
of the symmetric group Sn generate the whole group approaches 3/4 as n oo. 
Diaconis and Shahshahani [16] studied random walks on the symmetric group gen¬ 
erated by random transpositions and developed methods for estimating the rate of 
convergence. 

Since a folding map agrees with the identity map on one half-space, and with a 
reflection on its complement, it is a piecewise isometry. The dynamics of piecewise 
isometries have been studied in various contexts by numerous authors, for example 
see the survey by Goetz [17]. Different from interval exchange maps, our folding 
maps are neither one-to-one nor onto. 

On an abstract level, our results are motivated by classical theorems about di¬ 
rected graphs, and more generally Markov chains on discrete state spaces. A graph 
is transitive if every vertex x can be connected to every other vertex y by a path 
of directed edges. Finding a minimal subset of edges that still connects all vertices 
is known to be a hard problem, particularly if the graph has cycles. As a practical 
alternative, Aho, Garey, and Ullmann developed the theory of transitive reduc¬ 
tions, which are graphs on the same vertex set with the smallest possible number of 
edges [18]. In our setting, the points x S play the role of vertices, and a pair 
(x, u) with u £ G plays the role of a directed edge connecting x to F^x. Edge paths 
correspond to trajectories, and transitivity means that all trajectories are dense. 
Theorem 1 gives necessary and sufficient conditions for transitivity, Proposition 4.3 
constructs a transitive reduction, and Proposition 5.4 establishes the presence of 
cycles. The random walk defined in Eq. (2) corresponds to the Markov chain de¬ 
fined by a weighted directed graph, with ^ playing the role of the edge weights. 
Theorem 2 and Gorollary 2 concern the steady-state of the Markov chain. 

Another parallel can be drawn between Eq. (2) and shift dynamics. Suppose that 
the probability density ^ is concentrated on a finite set G of directions. Eor a given 
initial point, x £ if we omit from the underlying sequence of directions those 

terms where the folding map acts as the identity (rather than by reflection), we 
obtain a random sequence in G where only certain transitions are admissible. The 
conditions for admissibility, however, depend on the current location of the random 
walk in for this reason the sequence of directions is not simply a subshift of 

finite type. 


4. Conditions (Cl) and (C2) 

We now discuss the hypotheses of the main theorems. The geometric condition 
(Cl) can be expressed in a number of different forms. Note that only parts (c) and 
(d) of Proposition 4.1 are needed subsequently. 

Proposition 4.1. Let G C . The following are equivalent: 

(a) The open half-spaces Hu with u £ G cover 

(b) G is not contained in any closed hemisphere ofS‘^~^; 

(c) the convex hull of G in contains the origin in its interior. 

Let G be the support of a Borel probability measure p,. Then the above conditions 
hold if and only if 

(d) Q < < 1 for all x £ 
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Proof, (a) =?> (b): If G satisfies (a), then for every x S there is a direction 

u € G such that —x € Hu- This means that x ■ u < 0, implying that G is not 
contained in the closed hemisphere {u G | a; • m > 0}. Since x was arbitrary, 
this shows (b). 

(6) ^ (c): If (c) does not hold, then there is a hyperplane through the origin 
that does not meet the interior of the convex hull of G. Therefore, G is contained 
in a closed half-space, which intersects in a closed hemisphere, contradicting 

(b). 

(c) (a): Let G be the convex hull of G in Given a point a; in consider 

the linear functional defined by ^{y) = x ■ y on If the origin is an interior point 
of G, then £ takes both positive and negative values on G. Since £ assumes its 
maximum at an extreme point of G, and the extreme points of G are contained in 
the closure of G, we conclude that £{u) > 0 for at least one u G G. Thus x G Hu, 
establishing (a). 

(a) (d): Let ^ be a Borel probability measure on and x G Assum¬ 

ing (a), the support of y, contains a direction u with x G Hu, that is, u G H^. Since 
Hx is open, it follows that y{Hx) > 0, proving the first inequality in (d). Replacing 
x with —X shows that 1 — y{Hx) > y{H_x) > 0, proving the second inequality in 

(d). The converse implication is obvious. □ 

Denote by 

(G) := {Run ■ ■ ■ Rui S 0{d) I n > 0, ui,..., Mn G G} 

the subgroup generated by the reflections i?„ at hyperplanes with unit normals 
u G G. The algebraic condition (C2) can fail in three ways: If (G) lies in a lower¬ 
dimensional subgroup, if it splits into two subgroups that act on orthogonal sub¬ 
spaces, or if (G) is a hnite Coxeter subgroup of 0{d). Thus, assumptions (1) and 

(2) in the following proposition are also necessary. In dimension d > 2, assumption 

(3) is stronger than necessary, because only integer multiples of j, and ^ can 
appear as angles between elements of a finite Coxeter subgroup that acts on R'^ [19, 
Theorem 9 and Proof of Lemma 4.2]. For (d-l-l)-element sets G C a precise con¬ 
dition for generating finite Coxeter groups was obtained by Felikson [20]; for more 
general sets of reflections in R*^ this is an open problem. A set of conditions slightly 
stronger than (I), (2), and (3) was used by Eggleston [21, proof of Theorem 46] 
and Klain [22, after Corollary 5.4] to establish convergence to balls for sequences of 
Steiner symmetrizations. 

Proposition 4.2. Let G C If (1) G spans R'^, (2) G is not a union of two 

non-empty mutually orthogonal subsets, and (3) not all angles between directions 
are commensurable with tt, then the subgroup (G) is dense in 0{d). 

Proof. We argue by induction over the dimension. For d = 1 there is nothing to 
show. For d = 2, by (1) and (3) there are directions u,v G G such that the angle 
d(u, v) is incommensurable with tt. Since the composition RuRv generates a dense 
subgroup of rotations in SO{2), it follows that ({m,i;}) is dense in 0(2). 

Suppose now that the proposition holds in dimension d — 1, where d > 2. If 
G C satisfies (1), (2), and (3), then it contains a subset G' that spans some 
hyperplane in R'^ and also satisfies (2) and (3). Let Sy be the orthogonal group 
on z;-*-. By the inductive hypothesis, (G') is dense in Sy. Note that Sy has a unique 
pair of fixed points at iw, and acts transitively on the unit sphere in v^. 
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By (1) and (2) there is a direction u G G that is linearly independent but not 
orthogonal to G'. In particular, w = RuV is linearly independent of v. Therefore 
the conjugate subgroup Sw = RuSyR^^ is different from Sy Intersecting Sy and 
Sw with SO{d), we obtain two distinct subgroups conjugate to SO(d—l) x {1} 
in (G). But SO{d) has no non-trivial compact subgroup that properly contains a 
copy of SO{d—l) [23, Lemma 4]. It follows that (G) contains SO{d). Since it also 
contains orientation-reversing elements, (G) = 0{d). □ 

Our next results concern finite subsets of 

Proposition 4.3. The minimal number of directions inG G required to satisfy 
(Cl) and (C2) is d+1. 

Proof. Condition (Cl) requires at least d+1 points, the minimal number of extreme 
points for a convex set with interior in 

To construct a set of d -I- 1 points that satisfies both (Cl) and (C2), we start 
from a basis of unit vectors ui,..., for The convex cone generated by this 
basis has non-empty interior. We choose a unit vector Ud+i whose antipode —Ud+i 
lies in the interior of the cone, and, moreover, Ud+i encloses an angle with ui that 
is not a rational multiple of tt, and set G = {ui,..., Ud+i}- Then (Cl) holds by 
Proposition 4.1(c), and (C2) by Proposition 4.2. □ 

We do not have a good characterization of minimal subsets (under inclusion) 
satisfying (Cl) and (C2). Such minimal subsets may contain more than d -|- 1 di¬ 
rections. For example, both conditions hold for G = {ui, —ui,..., Ud, —Ud}, where 
ui and U 2 enclose an angle incommensurable with tt and ui,... ,Ud form a basis of 
unit vectors in but every proper subset is contained in a closed hemisphere and 
thus fails (Cl). 

There are sets in satisfying (Cl) and (C2) that have no minimal subsets. 
For example, the set G = | fc > 1} C satisfies (Cl) and (C2), and so 

does every infinite even subset, while every finite subset generates a finite dihedral 
subgroup of 0(1). Our next result shows that in higher dimensions, minimal subsets 
always exist, and are finite. 

Proposition 4.4. Let d > 2. If G C satisfies (Cl) and (C2), then there is a 
finite subset G' G G that also satisfies both conditions. 

Proof. If G is finite, there is nothing to show. If G is infinite, we choose a small 
d-dimensional simplex inside the convex hull of G in R'^ that contains the origin in 
its interior. By the Krein-Milman theorem, each vertex of the simplex is a convex 
combination of a finite subset of G. The convex hull of the union of these subsets 
contains the origin in its interior. By Proposition 4.1(c), it satisfies condition (Cl). 
To obtain G', we include up to two more elements of G in the union in such a 
way that G' does not split into mutually orthogonal subsets, and that it contains 
a pair of vectors enclosing an angle that is not an integer multiple of j, or 
In dimension d > 2, this implies that G' is not contained in a finite Coxeter group. 
Thus (Cl) and (C2) both hold for G'. □ 

5. Invariant subsets and dense trajectories 

This section is dedicated to the proof of Theorem 1 and Corollary 1. Let 
and Ry be the folding map and the reflection defined by a direction u G By 
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definition, a set A C is positively invariant under if 

FuiA) C A; 

if RuA = A then A is invariant under i?„. We say that A is almost positively 
invariant if FuA \ A is a null set with respect to the uniform probability measure 
(t; if the symmetric difference A A RuA is a null set then A is almost invariant. 

We consider sets that are positively invariant under the folding maps indexed by 
a non-empty set of directions G C Under condition (Cl), positively invariant 

subsets are almost invariant. 


Lemma 5.1. Let G C be a subset that is not contained in any closed hemi¬ 

sphere. If A C is almost positively invariant under F^ for all u € G, then it 
is almost invariant under i?„ for all u € G. 


Proof. If A is almost positively invariant under F^, then {RuA n Hu) \ A is a null 
set. In that case, 


/ x-uda{x)= / x-uda{x)— / x-uda{x)>0. 

Ja J ArH„ J{RuA)nHu 


Equality holds only if {AriHu)\RuA is a null set, in which case A is almost invariant 
under i?„. 

Fix Mo G G. By Proposition 4.1, the origin is an interior point of the convex hull 
of G. This means that for £ > 0 sufficiently small, —euo can be represented as a 
convex combination 


n 


(4) -euo = '^ aiUi 

i=l 

with Ui € G and > 0 for 1 < i < n. It follows that 


/ X ■ Uodx cti X ■ Ui da{x) = 0 . 

Ja Ja 


By the positive invariance of A, all summands are nonnegative. Since £ > 0, the 
integral f x ■ Uq da{x) vanishes, and therefore A is almost invariant under Rug- Cl 


Under condition (C2), the sphere has no non-trivial almost invariant subsets. 

Lemma 5.2. Let G C If (G) is dense in 0{d), then the only subsets ofS‘^~^ 

that are almost invariant under the reflections {Ru \ u € G} are null sets and their 
complements. 

Proof. The claim can be proved directly, but we give a shorter proof which uses 
spherical harmonics. Assume that A C is almost invariant under Ru for all 
u € G. Consider the indicator function of A as an element of L^{S‘^~^), and expand 
it in spherical harmonics as 

OO 

1 ^ 1 ( 2 ;) = '^Yk{x). 
k=0 

Here, Yk is orthogonal projection of 1 a onto the spherical harmonics of degree 
fc, and the series converges in the mean-square sense. Since this representation is 
unique, each Yk is invariant under composition with Ru for all u € G. Therefore, 
Yk is invariant under (G), and, by continuity, under the entire group 0{d). But 
for k > 0, the action of 0{d) on the spherical harmonics of degree k fixes only the 
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zero polynomial. It follows that Yfe = 0 for all fc > 0, and the constant function 
Yq = (j{A) agrees cr-almost everywhere with 1^- □ 

Combining Lemmas 5.1 and 5.2, we conclude that there are no non-trivial compact 
positively invariant subsets. 

Lemma 5.3. Assume G C satisfies (Cl) and (C2). Then no compact subset 
(and no open subset) other than 0 and is positively invariant under 

all foldings | u G G}. 

Proof. Let A C be a non-empty compact set that is positively invariant under 
Fu for all u € G. By Lemma 5.1, A is almost invariant under for all u € G, and 
by Lemma 5.2 either A or its complement is a null set. We want to exclude the first 
alternative. 

For every e > 0, 

A^ = {x & I d{x,A) < e} 

is a compact set of positive measure. Since A is positively invariant and the maps 
are non-expansive, A^ is positively invariant as well. By Lemmas 5.1 and Lemma 5.2 
its complement in is a null set; in particular, A^ is dense in By com¬ 

pactness, A^ = and hence A = f]As= 

If, on the other hand, A is a non-empty open set that is positively invariant under 
G, then its complement is a compact set that is positively invariant under —G. By 
the first part of the proof, it is empty. □ 

Proof of Theorem 1. For x G , consider the orbit 

G*a; = {Fu„ ■.. n > 0, ui,... u„ G G} . 

By definition, G*x contains all trajectories of x under sequences (Fu^) with direc¬ 
tions Un G G. Since G*a; is positively invariant under F^ for u € G, its topological 
closure is a positively invariant non-empty compact subset of By Lemma 5.3, 
the only such subset is the entire sphere. It follows that Gi,x is dense in Hence 

there exists for every e > 0 and every x,y G a finite sequence of directions 

ui,... Uji in G such that d{Fu„ ■ ■ ■ Fu^x, y) < e. 

We claim that the sequence can be chosen independently of x and y. In fact, for 
every e > 0 there is a finite sequence ui,.. .un in G such that 

(5) min d{Fu„ ■■.Fu^x,y) <e Vx,y . 

n<N 

The sequence is constructed by concatenating a finite number of shorter segments 

51., Sk as follows. 

Given £ > 0, cover by finitely many open balls Hi,..., Bk of radius e/3 
centered at ci,..., cx. To construct Si, choose a finite sequence of directions in G 
such that the corresponding trajectory starting at ci visits the ball B 2 , and then 
extend that sequence so that the trajectory visits each of the balls Bi,... Bk- The 
segments Sk for 1 < k < K are constructed inductively. Assuming Si,... ,Sk-i 
have already been chosen, let yk be the final point of the trajectory of Ck under 

51.. .., Sk-i- Choose Sk such that the trajectory of yk under Sk visits each of the 
balls Hi,..., Bk- Then the trajectory of Ck under Si,... ,Sk visits each of the balls. 

Let ui,... ,un he the sequence of directions given by 5i,..., Sk- For x,y € 
let Bi and Bj be the balls containing x and y, respectively. By construction, 

F F r G B 
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for some n < N. Since foldings are non-expansive, the triangle inequality implies 

d{Fu^ ■ ..Fu^x,y) < d{x,Ci) + . ..Fu^a.Cj) + d{cj,y) < e. 

This establishes Eq. (5). The desired infinite sequence (m„) is obtained by concate¬ 
nating the finite sequences constructed above for e = 2~^ with j > 1- □ 

Proof of Corollary 1. Let 0 be a continuous function on such that Su(j> = (j) for 
all u € G. Then cfoFu > (j> for all u G G, that is, (j) increases along trajectories. We 
need to show that the restriction of <j) to each centered sphere {|x| = i?} is constant. 
By scaling, it suffices to consider the case R= 1. 

Let x,y G 8"^“^ be given. By Theorem 1 there exists an infinite sequence of 
directions (u„) such that the trajectory (x„) defined by Eq. (1) is dense in 
Choose a subsequence {Xn^) that converges to y. By monotonicity and continuity, 

(j){x) < lim = (j){y). 

K—¥00 

Switching the role of x and y yields the reverse inequality <f>{y) < 4>{x). We conclude 
that (j) is constant on □ 

Theorem 1 says that there exists a sequence of directions in a set G (satisfying 
(Cl) and (C2)) that generates dense trajectories for all starting points. However, 
given G, it is not obvious how to explicitly find a dense trajectory. Different from 
reflections, periodic sequences of foldings never generate dense trajectories. 

Proposition 5.4. Let {un)n>i be « periodic sequence of directions, with Un+p = 
for some integer p and all n > 1. Then Eq. (1) has no dense trajectories, and at 
least one trajectory is periodic. If, in addition, G = {ui,... ,Up} satisfies (Cl) and 
(C2), then the periodic trajectory is non-constant. 

Proof. Consider the composition F := Fu^ o ■ ■ ■ o F^ as a map from Hu^ to itself. 
Since F is continuous and Hu^ is homeomorphic to a ball in by Brouwer’s 

fixed point theorem there exists a point x G Huj, with F(x) = x. By construction, 
the trajectory (Xn)n>o of x is periodic, and in particular not dense. If a;(j is another 
trajectory, then the sequence d(xn,x'^) is non-increasing, because folding maps are 
non-expansive. Let r = Yim.d{xmx'.^)- If r = 0, then the limit points of (a;(j) are 
precisely xi,... ,Xp. Otherwise, let y be a limit point of {x'f^p)k>i, and let (?/„) be 
the trajectory of y. By construction, each point on the trajectory of y is a limit 
point of the subsequence (a;(j^;.p)fc>i. Therefore, d{xn,yn) = r for all n > 0, that 
is, (y„) lies in the union of the d — 2-diinensional subspheres of radius r centered at 
the points xi, ... ,Xp. Since this is a null set, neither (y„) nor (a;(j) is dense in 

For the last claim, suppose that (a;„) is constant. Then the singleton {a;} is 
positively invariant under G. By Lemma 5.3, either (Cl) or (C2) fails. □ 

On the other hand, most trajectories are dense under the assumptions of Theo¬ 
rem 1. 

Proposition 5.5 (Random walks are dense). Let (Un) be an i.i.d. sequence of 
random direetions on with distribution y. If the support of y satisfies (Cl) 

and (C2) of Theorem 1, then almost surely the trajectory in Eq. (2) is dense in 
for every starting point x. 

Proof. Let s: > 0 be given. For n > 1, let be the event that for every pair of 
points x,y G there exists an integer n such that the random trajectory starting 

at X intersects an open e-neighborhood of y within some segment Xn+i, ..., X^+n- 
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By Theorem 1, there exists an TV < oo and a sequence (it„) in the support of ^ such 
that occurs on Xi, ... ,Xn- Since finite segments of trajectories depend con¬ 
tinuously on the sequence of directions, the probability of is strictly positive. 
By the Borel-Cantelli lemma, almost surely the event occurs infinitely often. 
Since e was arbitrary, the trajectory is almost surely dense. □ 

6. Random walks and invariant measures 

In this section, we prove the ergodicity results in Theorem 2 and Corollary 2. 
For the construction of the invariant measure, we need some more notation. 

Consider for the moment a single random direction U in with probability 
distribution /i. The random folding Fjj pulls a Borel function (p on back to 

(6) T^(j){x) = E{(j){Fu{x)) = ( (j){Fu{x))dn{u). 

The operator is linked to the action on measures defined in Eq. (3) by the 
change-of-variables formula 

[ {Tf,(p)dv= [ (j)d{Tf,#iy). 

Clearly, is a positivity-preserving linear operator that fixes constant functions. 
Since folding maps are non-expansive, preserves or improves the modulus of 
continuity of a continuous function. 

The following lemma establishes a useful monotonicity property for the extrema 
of (p and Tfj^p. 

Lemma 6.1. Let ^ be a probability measure on 8"^“^ whose support satisfies (Cl). 
Define by Eq. (6), and let p be a continuous function on 8"^“^ with maxp = M. 
Then maxT^c/) < M, and for any x G 

T^Pix) =M p{x) = M . 

The corresponding statements hold for the minimum of p. 

Proof. Left x G Since Fu{x) = x for all u G Hx, and p{x) < M for all u, we 

have 

Tfj_p{x) = / p{x) dp{u) + / p{Rux) dp,{u) 

< P{x)pl[Hx) + Mp{E,‘^~^\Hx) 

< M. 

Since p,{Hx) > 0 by Proposition 4.1, equality implies p{x) = XI. □ 

The random walk is associated with a canonical Markov process on 17 = nr!,>o 
the space of sequences on endowed with the product topology. If $ 

is a function on 17, we write 

Exm{X„)n>o)) = J {iFu„ ■ . ■ Fu,x)n>o) dp® 

for the expected value of 4) on the random walk starting at Xq = x. Here, 
p® is the product measure that defines the distribution of the sequence (t/„). Cor¬ 
respondingly, the probability of an event H C 17 is denoted by Px{A) = Ex{1a)- 
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The canonical Markov process is completely determined by either of the opera¬ 
tors Tfj, or By the Markov property, 

The next result says that converges uniformly to a constant function. 

Proposition 6.2. Let ^ be a measure on whose support satisfies condi¬ 

tions (Cl) and (C2) of Theorem 1. For every continuous function (p on 
there exists a constant <f> such that 

lim Ex {(piXn)) = ^ 

n—^oo 

uniformly in x. Moreover, 

min 4> < 4> < max (p . 

Both inequalities are strict unless p is constant. 

Proof. Let || - || denote the sup-norm on the space of continuous functions. Set 
pn{x) = Ex{p{Xn)). Since pn = (T/x)"')', its norm ||(()„|| decreases monotonically 
and the modulus of continuity of pn improves with n. By the Arzela-Ascoli theorem, 
there exists a subsequence {pnk)k>i that converges uniformly to some limiting func¬ 
tion, p. After passing to a further subsequence, we may assume that the sequence 
of gaps Uk — Uk-i increases strictly with k. 

We want to show that p is constant. Consider the sequence {'pn) defined by 
pn = Clearly, 

Wpm - Pm-\-n\\ = ||(T)x)’"(<^ “ Pn)\\ < \\P - Pn\\ 

for all m, n > 0. We use the triangle inequality and set m = Uk — rife-i, n = Uk to 
obtain the bound 

WPnic-nk-i ~ ^11 ^ WPnk-nk-i ~ Pnt 11 + I \Pnk ~ P\\ 

<11^ — pnf,_i 11 + I \ Pnk ~ ^11 j 

which converges to zero by the choice of the subsequence (nfe). It follows that the 
subsequence pn^-uk-i converges uniformly to p. 

Let M = max((). Since m.&'x.pn is non-increasing in n, and a subsequence con¬ 
verges to pa = p, we must have maxpn = M for all n > 0. By Lemma 6.1, the sets 
An = {x : pn{x) = M} form a decreasing chain. Their intersection A = p| A„ is a 
non-empty compact set that is positively invariant under E^ for /i-a.e. u S 
By Lemma 5.2, A = But this says that pn = M for all n, and the same holds 

for their limit p. Since the sequence \ \pn ~ P\\ is non-increasing, the convergence of 
the full sequence follows from the convergence of the subsequence. 

Clearly, min p < p < max p, since min pn is non-decreasing and max pn is non¬ 
increasing. Let M = max^, and consider the decreasing chain of compact subsets 
An = {x : pn{x) = XI}. Since the intersection A = f]An is compact and positively 
invariant, by Lemma 5.3, it is either equal to 8"^“^ or empty. In the first case, 
p = M is constant. In the second case, by compactness, A„ is empty and thus 
Tiiaxpn < M for some sufficiently large n. By monotonicity, p < maxpn < M. The 
same argument shows that min^ < p unless p is constant. □ 

Proof of Theorem 2. Let p he a continuous function on and let (W„) be the 

random walk defined by Eq. (2). By Proposition 6.2, Ex {p{Xn)) converge to a 
constant function, p, uniformly in x. It is a fact of Ergodic Theory that this 
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uniform convergence is equivalent to the unique ergodicity of the Markov chain [24, 
Theorem 4.10]. We show the part of the proof that we need here. 

The map ^ ha (/> is linear and continuous with respect to the topology of uniform 
convergence, and its value on nonnegative functions is nonnegative. By the Riesz- 
Markov theorem, there is a unique regular Borel measure p such that 

^ = / (jidp. 

Jsd-l 

Since the constant function (^ = 1 is mapped to itself, p is a probability measure. 
We next verify that p is invariant. For every continuous function (j) on 


(j>d{T^#p)= {T^(j)) dp = . 


Since 

T^4>= lim Ex{Tf^(j)[Xn)) = lim Ex{4>{Xn+i) = (f) = [ (t^dp, 

n—¥cci n—¥oo 

the measure represents the same distribution as p. By uniqueness, = p. 

It remains to show that the support of p is Given an arbitrary non-empty 

open set A C let be a nonnegative continuous function supported on A that 

takes values in [0,1] and does not vanish identically. Then p{A) > ^ > 0 by the last 
part of Proposition 6.2. Since A was arbitrary, the proof is complete. □ 


As an immediate consequence of the uniform convergence proved in Proposi¬ 
tion 6.2, we obtain the mixing property 

lim [ Ex{(i){Xn))ip[x)dp{x) = 
n—^oc Jgd-i 

for every pair of continuous functions (f),ip on 


Proof of Corollary 2. Let p and p be as in Theorem 2. We will show that for every 
p-integrable function on 

1 f 

(7) lim (t>[y)dp{y) 

almost surely for p-almost every x € The claim then follows by setting 

(f = 1 a - 

The proof of Eq. (7) calls for a standard application of Birkhoff’s ergodic the¬ 
orem to the canonical Markov chain associated with the random walk. See for 
example [25, Chapter 6]. The invariant measure induces a probability measure p* 
on n by 

p*{A)= f P4A)dpix). 

This measure is invariant under the left shift L{{Xn)n>o) = (^n+i)n>o, and every 
shift-invariant subset A C has p*(A) = 0 or p*{A) = 1. By Birkhoff’s theorem. 
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for every p*-irLtegrable function $ on and for p*-almost every sequence {Xn)n>o- 
In the case where <&((X„)„>o) = depends only on the initial point, we have 

$(L'=(X„)„>o) = , [ ^dp*= [ ^dp, 

which yields Eq. (7) except for sequences (Xn) in a set _B C of /5*-measure zero. By 
definition of p* and Fubini’s theorem, Px{B) = 0 for p-almost every x € □ 

If (/) is a continuous function on then the functions 

uniformly equicontinuous in x for all > 1 and every sequence of directions (C/n), 
see Eq. (2). Since is separable, it follows that Eq. (7) almost surely holds for 
every x G 


7. Properties of the invariant measure 


Finally, we study the properties of p. We find that p is generally not the uniform 
measure on the sphere. 

Proposition 7.1. Under the assumptions of Theorem 2, the invariant measure is 
uniform on if and only if p is even under u i-A —u. 

Proof. Recall that tr denotes the uniform probability measure on We want to 

show that 

T^#cr = a <;=> p{A) = h'i—rU) for all Borel sets A C . 

=>: Suppose Tf_i#a = a. Then, for every v G and e > 0, 


[ cr{Fu^iBe{v)))dp{u) = a{B^) . 


/S<i- 

By definition of the folding map, the integrand is given by 
a{F-\B,{v))) = 2a{B,{v)nHu) . 
We divide by cr^B,;) and take e —>■ 0, 

2a(R,(u)ni7„) [ J 

lim -;-^- = iim -;-7- = < i 

£->■0+ <^{Be) e^0+ <^{Be) 1 q 


if u • u > 0 , 
if u • u = 0 , 
otherwise. 


Integrating both sides over Hy and using that = cr, we obtain by dominated 

convergence 

a{Fff\B,{v))) 


1 = lim 
£->■ 0 + . 


■ dp{u) = 2p{Hy) + p{dHy). 


1 Cr(Be) 

It follows that p{Hy) = p{H-y) for all v G By Lemma 7.2, which is proved 

below, p is even. 

<^=: Conversely, if p{A) = p{—A) for all A C then 
T^#<j(A) = I aiF-\A))dpiu) 

= \j <r{F-\A))+a{FZ:,{A))dp{u) 

-if 


a'(A) + a{RuA) dp{u) 
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In the second line, we have used that ^ is even to change the variable u to —u in half 
of the integral. The third line follows, since F^^{A) contains a copy of A fl H±u 
together with its mirror image. In the last step we have exploited the reflection 
invariance of cr. □ 


Lemma 7.2. Let ^ be a regular Borel measure on If 

KHu) = KH-u) 

for all hemispheres Hu with u S then pL is even. 

Proof. Assume for the moment that p, is absolutely continuous with respect to the 
uniform measure a on with a smooth density m{x), and consider the expansion 
in spherical harmonics 

m{x) = '^Yk{x). 

k>o 

By the Funk-Hecke formula [26, Theorem 9.7.1], there exist constants such that 


( 8 ) 


Yk{x) da{x) = jkYkiu) 


'Hu 


for all M S In particular. 


p{Hu) = / m{x) da{x) = ^ 'jkYkiu) ■ 

k>0 


Since p{Hu) = p{H-u) for all u G and Yk{u) = {—l)^Yk{—u), the component 
yfcYfc must vanish for each odd k by the uniqueness of the expansion. We claim that 
7 fc 0 for odd fc, forcing Yk = 0. 

A key point of the Funk-Hecke formula is that the constant yt depends only on 
the order (fc) of the harmonic polynomial. It can be computed from Eq. (8) by 
setting u = Cd, and choosing Yk to be the normalized zonal harmonic Zk, whose 
restriction to depends only on the last variable Xd- For fc > 0, we use that Zk 
is an eigenfunction of the spherical Laplacian with an eigenvalue — Afc < 0 to obtain 

f Zkda = ( AZkda = Xf^ ( -dx^Zkda. 

J{xd>0} J{xd>0} J{xd=0} 

In the second step, we have applied the divergence theorem on the sphere. Since 
Zk depends only on the variable Xd, the normal derivative that appears in the last 
integral is constant on the equatorial sphere {xd = 0}. 

For k odd, Zk vanishes on {xd = 0}. Since the eigenvalue-eigenvector equation 
for Zk is a homogeneous linear second-order ordinary differential equation in Xd^ 
the normal derivative cannot vanish simultaneously, and so the integral is non-zero. 
It follows that 

7fc = ,7 / N [ Zkda {k odd). 

Zk[&d) Jfx 

We conclude that 


'{a;d>0} 


p{A) = ^ f Yk da = pi-A) 

k>0 even ^ 

for all Borel sets A C proving the claim when p has a smooth density. 

Otherwise, we approximate it with smooth measures pe, defined by 


d'eiA) = [ piQAifjeiQ) daiQ), 
Jso{d) 
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where tps is a smooth probability density supported on an £-neighborhood of the 
identity in SO{d), and a is the uniform measure. Let u G Since Q{Hu) = 

Hqu = —QH-u by linearity, the assumption on /r implies that 

= KQH-u) {Q e SO{d)). 

Therefore, for all u G By the first part of the proof, fig is 

even. Taking e —>• 0 we see that fi is even as well. □ 

We conclude with some open problems. 

Question 7.3. Given a measure /r on whose support satisfies (Cl) and (C2), 
let p be the invariant measure from Theorem 2. Is it true that 

p{A) =0 a{A) = 0 ? 

We suspect that the answer is affirmative. To motivate this, we write Eq. (3) as 

Tij_#v{A) = [ [ 1 a{Fu{x)) dp,{u)diy{x), 

dS'i-i JS'i-i 

and observe that is absolutely continuous with respect to cr, whenever either 

p or is absolutely continuous. In particular, the absolutely continuous component 
of the invariant measure p is itself invariant. By the uniqueness part of Theorem 2, p 
is either absolutely continuous or purely singular. If p has an absolutely continuous 
component in its Lebesgue decomposition, then p is absolutely continuous. If p is 
singular, nothing is known. If p were singular as well, then the trajectories of the 
random walk would accumulate on a set of measure zero. 

Another problem is to characterize how quickly the state of the Markov chain 
associated with the random walk in Eq. (2) converges to the steady-state. 

Question 7.4. Given a measure p on 8"^“^ whose support satisfies (Cl) and (C2), 
fix X G and define the random walk Xn by Eq. (2). At what asymptotic rate 

does the distribution of Xn converge to p as n ^ oo? Is there a cut-off phenomenon? 

One approach is to analyze the convergence of fn = to its limit f = J f dp 

for continuous functions f on In the case where p = a is the uniform measure 

on the sphere, computations similar to those in [1, Proposition 5.2] and [27] suggest 
that the difference Wfn — f\\ decreases with n~^ in the number of steps. Can this 
rate of convergence be improved by a judicious choice of p, perhaps supported on 
a finite set? 
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